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Single-file diffusion refers to an infinite system of interacting particles diffusing on a 
line, or hopping on a one-dimensional lattice, which cannot overtake each other. Such 
crowded motion arises in many quasi one-dimensional systems (see e.g. 0il)- The 
name ‘single-hle’ is particularly apt if one studies the evolution of a tagged particle. 
(The tagged particle is indistinguishable from the host particles, it merely carries a 
mark.) The behavior of the tagged particle is drastically different from Brownian 
motion [7 14 . This has also been demonstrated is many experimental realizations 
of single-hle diffusion 15 17 . The goal of this article is to highlight these differences 


by comparing dynamical properties of the tagged particle in single-hle dihusion (known 
and those we compute below) with corresponding features of Brownian motion. 

We limit ourselves to the simplest version of single-hle dihusion with uniform initial 
density. Denote by X{t) the displacement of the tagged particle. There is no drift in the 
case when the density prohle is uniform, so the average displacement vanishes, {X) = 0, 
if we set W(0) = 0. Higher odd moments also vanish: (X^) = (X^) = ... =0. The 
second moment exhibits a sub-dihusive growth 


{X^) = 2V{p) Vi 


( 1 ) 


for large t. This remarkable behavior was originally established 50 years ago for 
the simplest single-hle system composed of Brownian particles with hard-core repulsion. 
Intriguingly, the same sub-dihusive growth ([^ applies to all single-hle systems, each 
one is parametrized by a single amplitude D{p) depending on the (uniform) particle 
density p. This amplitude can be expressed through the two fundamental transport 
coefficients, the dihusion coefficient D{p) and mobility cr(p), underlying the large-scale 
hydrodynamic behavior of single-hle dihusion: 

^(P) 

P^\/VkD{p) 

For a random process without drift, the second moment does not generally 
characterize its probability distribution. In the case of single-hle dihusion, however, 
Arratia proved that it does and he showed that the probability distribution of the 
tagged particle position converges to a Gaussian [^. This striking result has a strong 
psychological inhuence, one starts to think that the behaviors of the tagged particle and 
Brownian motion are qualitatively similar and after replacing t in formulas for Brownian 
motion by Vi one would recover correct results for single-hle dihusion. This intuition is 
wrong as we shall demonstrate by comparing a few basic features of Brownian motion 
with their single-hle analogs. 

The crucial general distinction between Brownian motion and single-hle dihusion is 
that for the latter one must carefully dehne the averaging. The dynamics is stochastic, 
so we obviously perform the averaging over diherent realizations of dynamics. Another 
source of randomness is hidden in the initial conditions. Single-hle dihusion is an inhnite 
particle system, so even if the tagged particle is initially at the origin, X(0) =0 as we 
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Figure 1. Plot of R{r) in Q for the two initial settings. The function R{r) does not 
depend on T>{p) and increases monotonically from 0 to 1 with increasing r. 


always assume, the initial positions of other particles affect the behavior of the tagged 
particle. If we consider a fixed initial condition which is macroscopically uniform (e.g., 
the equidistant distribution), instead of Q one finds 10-12 

(X2)det = V{p) V2t (3) 

The usually cited sub-diffusive growth Q tacitly presumes that we additionally perform 
the averaging over all (macroscopically uniform) initial conditions. It is not surprising 
that the behaviors ([^ and (|^ are different for small and moderate times. In the 
long-time limit, however, one anticipates that an equilibrium distribution of particles 
is reached even if we start with a very non-equilibrium, e.g. equidistant, initial 
distribution. This seemingly implies that the initial condition will be ‘forgotten’ and the 
mean-square displacement (and any other characteristics, e.g., (X^^)) should not depend 
on the averaging over the initial conditions. What looks even more striking is that for 
any initial condition the mean-square displacement has the same asymptotic behavior, 
Eq. (^, yet after averaging we obtain the result which is y/2 times larger, Eq. ([^. 

n this paper we study a few basic properties of single-file diffusion that have simple 
counterparts in the context of Brownian motion. One such property is the two-time 
correlation function. The neat formula 


{x{ti)x{t2)) = 2D min(fi,f2) 


( 4 ) 


for the two-time correlation function encapsulates the major properties of Brownian 
motion, e.g., it can be derived from the equal-time correlation function {x"^) = 2Dt 
using Markov property. For single-hle diffusion, we shall obtain (for large ti and ^ 2 ) 


(X(fi)X(f2))ann 


V{p) 


y/ti -|- \/t2 — \/|t2 — ti 


( 5 ) 


if we perform the averaging over the initial conditions and 


{W(fi)W(f2))det = V{p) 


Vtl + ^2 — ■\/|^2 — 


( 6 ) 


for a fixed initial condition. The derivation of Eqs. (§-(§ for an arbitrary single-file 
system is one of our major results. To illustrate the difference between the two settings 
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we plot (figure the ratio 


(■Y(ti).Y(i;)> 



(7) 


This quantity is also known as the Pearson correlation coefficient. 

deduce ([^ from (@. Indeed, let us allow the evolution 


Intriguingly 18 


one can 


to proceed up to time T, and start measuring the displacement of the tagged particle 
after that moment. Denote by T + ti and T + ^2 two later times and by 


Y(tj) = X{T + t,)-X(T), ] = 1,2 (8) 

the proper displacements. The two-time correlation function is 
{Y{tfiY{t,)) = {X{T + tfiX{T + 

^2))det 

- {X(T + «i).Y(r))*. - {X(T + (2).Y(r)>de, 

Using ([^ we re-write {Y{ti)Y(t 2 )) as 

— 2T Y ti y/ti — \/ 2T -1- ^2 + 


We now take the 


to 

^ 0 , ^ 0 

T ’ T 


limit and hnd 


(9) 


IfY(ti)Y(t2)) — D{p) a/U -|- \/t2 ~ \/|t2 — ti 


( 10 ) 


Thus we do not need to average over initial conditions: Performing the shift of time 
(|^ and moving far ahead, as in Eq. (§. we effectively put the system into equilibrium, 
so should be the same as ([^ and indeed it is. Needless to say, ([^ and ([^ reduce 
to ([^ and ([^ for equal times and this explains why the amplitudes in (|^ and ([^ are 
different, and why the difference is universal by a factor ^/2. 

Many quantities arising in the context of Brownian motion are described by the 
arcsine law. This law gives, e.g. the distribution of time that the Brownian particle 
spends to the right of the origin and the distribution of time at which the Brownian 
particle reaches maximum. Our second goal is to investigate their single-hle analogs. 
Consider the trajectory {X(t) : 0 < t < T} of the tagged particle and denote by 


T+ = Measure({t G [0,T] ■ X(t) > 0}) 


( 11 ) 


the total time that the particle spends to the right of the origin. This time scales 
linearly with T, so one expects the cumulative distribution to be the function of the 
ratio T+/T, that is, Prob(T+ < rT) = Pi(r) in the T —)■ cx) limit. Similarly denote by 
T* the time at which the displacement of the tagged particle reaches its maximum. This 
time T* is a random quantity which scales linearly with T. Hence we anticipate that 
Prob(T* < tT) = P 2 ('r) in the T —)■ cx) limit. 
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For Brownian motion, the distributions Pi(t) and P 2 (t) coincide and are known as 
arcsine distributions (Toj-j^: 

2 

Pi(r) = P 2 (r) = — arcsin(A/r) (12) 

TT 

Finding the distributions Pi(r) and P 2 ('r) for single-hle diffusion is an interesting 
challenge. We have computed one non-trivial moment (r^) = dr The 

result is independent of the implementation of the single-file constraint, but depends 
on whether we perform or not the averaging over initial conditions, and both answers 
are smaller than the corresponding result (r^) = | for Brownian motion. 

The third characteristic we study is the growth of the maximal displacement. Recall 
that for Brownian motion the typical spread grows as y/lH, but maximal fluctuations 
grow slightly faster. It is useful to think about the smallest asymptotic envelope that 
contains Brownian trajectories. Informally, one can write 

\x{t)\ < \/2{x‘^{t)) In(lnt) = ADt In(lnt) (13) 


This is often referred to as the law of the iterated logarithm 19,^. The precise 
meaning of (13) is that |x(t)| < C ^/ADt\niJxlt) is asymptotically correct (i.e., valid for 


sufficiently large time) for any C* > 1, whereas for any C < 1 any Brownian trajectory 


crosses the boundaries of envelope for arbitrarily large times. More formally 19| 21| , 

xit) 


lim sup —-p _ 


= 1, lim inf , 


= -l. (14) 


The extension to the single-file diffusion is rather straightforward. We shall argue 


that the analog of (14) is 


Xit) 

hm sup — , 

v/(X2(t)) Inf 


= 1, hm inf — 


= - 1 . 


(15) 


The choice of (X'^it)) depends on whether we perform or not the averaging over initial 
conditions. 

Equal time and two-time correlation functions of the tagged particle displacement 
[Eqs. (0) (|^) foil us a lot about single-file diffusion. Accordingly, one would 

like to compute multi-time correlations, and more generally the probability of an entire 
trajectory. For Brownian motion, the latter probability can be represented by a path 


integral. In this formalism, the probability of a trajectory is expressed as 22,23 


P[x] ~ exp{—A[a;(r)]} 


( 16 ) 


where the functional of the Brownian trajectory 

2 






'0 


\dT 


(17) 


is known as action as it is formally identical to the action of a classical free particle. 

For single-file diffusion, there is no comparably simple formula for the action. A 
macroscopic fluctuation theory (MFT), a recently developed systematic approach to 
non-equilibrium fluctuations in diffusive systems (see 24-26 for a review), allows one 
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to proceed. Using MFT we have shown 11 how to analyze the full statistics of the 
tagged particle position in a general single-hle diffusion. Here we extend the formalism 
to assign probability weight to an entire trajectory of the tagged particle. In particular, 
we study a Legendre transform of the action and express it in terms of solution of a 
variational problem. At a formal level, the analysis boils down to solving Euler-Lagrange 
equations corresponding to the variational problem. A series expansion of the Legendre 
transform generates all the multi-time cumulants of tagged particle position. 

The advantage of our approach is its generality, viz. it applies to all single-file 
systems. The microscopic details of the system, e.g. the inter-particle interaction 
and the size of the particles, appear at the macroscopic scale only in terms of two 
transport coefficients D{p) and cr(p). It is still extremely difficult to solve the Euler- 
Lagrange equations for an arbitrary single-file system, essentially the only fully solvable 
case corresponds to a single-file system composed of Brownian particles. Less complete 
results, however, can be established for an arbitrary single-file system. Specifically, we 
shall describe a perturbation scheme which allows, in principle, to compute the multi¬ 
time correlations order by order. Using this method we derive general formulas (|^-(|^ 
for the two-time correlation function. 

One important result that emerges from our macroscopic analysis is that the multi¬ 
time cumulants of tagged particle position, knifi, ■ ■ ■ ,tn) = {Xifi) ■ ■ ■ X{tn))c) scale 
sub-diffusively. More precisely, all multi-time cumulants scale as 

knifi,- ■ ■ fin) ^ ^ (18) 

when T = max(fi,..., is large. Correlation functions {Xfii) ■ ■ ■ Xfin)) can be 
expressed in terms of the cumulants. Using these expressions one finds that for n > 2, 
the asymptotic behavior of the n-time correlation function is determined by the two- 
time cumulant (the one-time cumulant is zero due to symmetry), the contribution from 
other cumulants can be ignored. For example 

{X{ti)X{t2)d^itfi)d^itfi)) = ^2(U; ^2)^2(^3) ^ 4 ) + ^2(U; ^3)^2(^2; ^ 4 ) + ^2 (U ; ^ 4 ) ^2 (^2; ^s) 

+ kfitlfi2fi‘ififi) (19) 


and using (18) we see that the ^4 term is negligible compared to the rest in the T ^ 00 
limit. 

The fact that the asymptotic behavior is expressed through the two-time cumulants 
implies an important feature of the joint distribution of tracer position at multiple times, 
viz. this joint distribution is asymptotically a multivariate Gaussian distribution: 


P{X,,---,Xn)- 


1 






( 20 ) 


Here Xj = Xfij) and we used a vector notation X = (Xi, • • • ,X„) with X* being the 
transpose of X. Further, S is an n x n covariance matrix with = {X{ti)X{tj)) 
and |S| denotes determinant of S. (The covariance matrix depends on the setting, one 
should use ([^ and (|^ for the annealed and the deterministic cases, respectively.) The 
asymptotically Gaussian distribution of X{t) is a special case of (20). 
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A direct consequence of (20) is that all multi-time correlations can be expressed in 


terms of the two-time correlation using the Wick’s theorem (or the Isserlis’ theorem). 

The outline of this work is as follows. In section]^ we consider the time spent by the 
tagged particle to the right of the origin and compute the second and third moments. In 
section]^ we give an heuristic derivation of the ‘law of the simple logarithm’, Eq. (15). 
The MET framework is extended to the case of a functional of the entire trajectory in 
section and the formalism is used to derive the two-time correlation functions 
in section We summarize our findings in section 


2. The arcsine law 


Consider a Brownian motion x{t) which evolves during the time interval 0 < t < 1. Let 
T be the time spent by this x{t) on the positive side i.e., 

T = [ dtQ[x{t)] ( 21 ) 

Jo 

where 0[a:] is the Heaviside step function. The time r is a random variable depending 
on the history of the Brownian particle. The cumulative distribution P(r) is given by 
the celebrated arcsine law (12), see e.g. 19-21 . Equivalently, the arcsine distribution 
has the density 

= I rJ — V P2) 

^ V^(l - 

This formula shows that, contrary to a naive expectation, it is more likely that r is 


either small or close to 1, and it is least likely to be around r = 1/2. Using (12) or (22) 
one can compute all moments: 


^ 1 r(n+|)r(l) _ (2n-l)!! 

' ' 7r r(n + l) (2n)!! 


(23) 


The first moment, (r) = merely reflects the right-left symmetry; the following 
moments, e.g. (r^) = | and (r^) = cannot be deduced merely from the symmetry. 

Finding analogs of (12) or (22) is of course the main challenge. Here we present 
a detailed computation of the second moment of r. We then briefly explain how 
to compute the third moment. These calculations are pedestrian and laborious, and 
generalizing them to the forth moment already appears very challenging. In particular, 
we cannot compute (r”) for large n which would have shed the light on the structure of 


the single-hle analog of the density (22) 


Scale invariance of the Brownian motion allows one to study the behavior on an 
arbitrary time interval, so the common choice is 0 < f < 1. In the case of single-file 
diffusion we must consider long intervals: 0 < f < T with T 3> We shall 

measure all time variables in units of T, so e.g. r = T+/T is the fractional time that 
the tagged particle spends to the right of the origin. 

The average time which the tagged particle spends to the right of the origin is the 
same as the time it spends to the left, so (r) = 1/2. The first non-trivial moment is (t^). 
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This second moment is different from that of the Brownian motion. We will see that 
the result depends on whether we perform or not the averaging over initial conditions. 
The former setting is known as annealed, the latter is deterministic or quenched. 


Using the definition (21) we find 

(0[X(ti)]0[X(t2)]) 


(r^) = 2 r dh f 

Jo Jo 


(24) 


To complete the computation it is sufficient to know the joint probability of positions 
at two times. Let P(X, t) be this joint probability, where X = (Xi, X 2 ) and t = (^ 1 ,^ 2 )- 
From (20) we hnd that P(X, t) is Gaussian which can be written as 

1 


F(X,t) = 


27rcri(J2A/l — 


exp 


1 


where 


at = {Xt), ai = {X^), R = 


2(1 -P2) 

(X 1 X 2 ) 


Xf X1X2 xl 

-2R + -4 

af (Ti(T2 a 2 


(Ti(T 2 


(26) 


(26) 


These quantities are different in annealed and in deterministic settings: in the 
annealed case we have, using ( 1 ), ct| = D{p) when the initial condition is fixed we 

have, using (|^, cr| = D{p) \/^- The quantity R depends on the ratio r = ti/t 2 . From 
we get 


P(r) = 


1 + a/t - \/r 
2 rip 

a/1 + r - ^/Y 

y2rV4 


annealed 


deterministic 


(27) 


Using these results we can compute the average ( 0 [X(ti)] 0 [X(t 2 )]) in (24). One 
can avoid calculations by relying on a neat formula p. 937] 


(0[X(U)]0[X(t2)]) = / dX, 


1 1 

dX 2 P{X,t) = 7 + — arcsinP(r) 
4 zn 


(28) 


which is valid for any bivariate Gaussian distribution (25). Interestingly, the final result 
(28) does not depend on cxi and (J 2 , it depends only on the Pearson correlation coefficient 
R{r). For completeness, we note that Eq. (28) can be easily derived without referring 


to [^, it suffices to compute the Gaussian integral over X 2 and then perform the 
integration over Xi using an identity 


dx e~ 


'0 


erfc(gx) = arccos 

'71 


Q 


a /1 + 


Substituting (28) in (24) leads to 

/■i /-fe 


{r^) = 


dU 


dti 


'0 


'0 


- H— arcsinP(r) 

2 TT 


Recalling that r = ti/t 2 we simplify (29) to 

1 I 

(r^) = 7 H-/ dr arcsinP(r). 

4 27r ./n 


(29) 


( 30 ) 
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This together with (27) completes the calculation of the second moment. Numerically, 

0.353766 ... annealed 

0.326602 ... deterministic. 




(31) 


Equation (30) applies to any stochastic process for which the joint probability 


distribution P(X, t) is Gaussian. For instance, for standard one-dimensional Brownian 


motion R{r) = ^/r and plugging this into (30) we obtain (r^) = | which can be 


alternatively deduced from (22). 


It is remarkable that the final results for (r^) are universal, i.e. they do not depend 
on the details of the single-hle system. This is a general feature, it continues to hold for 
the higher moments as we show below. 


Is it feasible to use the same strategy for computing higher moments? Using (21) 
we obtain 

(r") = n! • • • r dt, (0[X(tO] ■ ■ ■ 0[X(f„)]) (32) 

Jo Jo 

so one should be able to compute the average (0[X(ti)] ■ ■ ■ 0[X(t„)]). This quantity 
is often referred to as the orthant probability [^. Some results about the orthant 
probabilities are known in the interesting to us situation when the joint probability 


distribution of Xitjfis is a multivariate Gaussian (20). From the general structure 


of (20) we already conclude that the orthant probabilities depend on the correlation 

cr?. Since 


coefficients Rij = R{ti/tj) and do not depend on (non-universal) variances 


j ■ 


R{r) are same for all single-file systems, all the moments (32) are universal. 

It is difficult, however, to establish explicit results despite of the fact that the joint 
probability distribution is a multivariate Gaussian distribution. Explicit (and simple) 


formulas have been derived only for n = 2, Eq. (28), and for n = 3. In this latter 


case 


(0[X(ti)]0[X(f2)]0[X(f3)]) = ^ [arcsinPi2 arcsinP23 + arcsinPis] (33) 

8 dvr 


Gombining (32) and (^33| we find 


r^) = —dr arcsinP(r) = 
8 dvr ./n 


0.280649... 
0.239903... 


annealed 

deterministic. 


(34) 


We also note that the second moment has been recently computed for the 
single-hle system composed of Brownian point particles in a very inhomogeneous setting, 
namely when all particles are initially at the origin. The result is (r^) = | + |-fo(7r) = 
0.33696..., where Jq is a Bessel function. The methods of Ref. 29 crucially rely on 


the composition of the system and hence they are entirely different from the methods 
which we use. 


3. Law of the simple logarithm 

Typical displacements in Brownian motion grow as \/i with time. Of course, there 
are excursions beyond the typical distance. A way to quantify the likelihood of these 
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excursions is to find the smallest envelope that contains the Brownian motion with 
probability one at infinite time limit. The answer is given in the law of the iterated 
logarithm, Eq. (14). Informally, the smallest envelop grows ■\/ln(lnf) times faster than 
the typical displacement. 

Here we argue that for the single-file diffusion the minimum envelope grows ■^ln(t) 
times faster than the typical displacement. Therefore the tagged particle has relatively 
wider excursions than the Brownian particle. 

We shall use heuristic arguments, so we first show that they lead to the well-known 
correct result for the Brownian motion. Let us think about the envelop as an absorbing 
boundary moving to the right according to 

b{t) = \f2a‘^{t) a{t) (35) 

where a{f) is a function that diverges as f —>■ cx) and cr^(t) = 2Dt is the variance. 
Denote by S{t) the survival probability of the Brownian particle. Since the envelope 
moves faster than the standard deviation, we assume that the probability distribution 
still has a Gaussian form, but with normalization S{t). We write 

S{t) 


P{x,t) = 


2a2(t)J 

The survival probability has to be determined self-consistently. This comes from the 
condition that the rate at which the S{t) decreases must be equal to the probability 
current at the position of absorbing envelope: 


X 


(36) 


dS _ 8P 

dt dx x=b(t) 

Plugging (36) into (37) we obtain 
dS _ _DS e 
dt 


a 


TT 


(37) 


(38) 


For Brownian motion = 2Dt. We insert this into (38) and note that if 

e~°‘y/a ~ (lnf)“^“'^ for any 5 > 0, then the survival probability approaches a finite 
value in the f — ?• cx) limit; however, if yW ~ (Inf)"^"'''^ for any <5 > 0, then the 
survival probability vanishes when f —)■ cx. The borderline case is thus obtained for 
^/a ~ (lnf)“^, i.e. a ~ In(lnf). This together with (35) leads to (14). 

For the single-file diffusion we proceed as in the case of the Brownian motion. 
At large times the distribution for the tagged particle position X(t) is asymptotically 
Gaussian, but with a variance that grows as cr^(f) ~ ^/i. Gombining this with 
(38) we hnd that the minimal growth rate of the envelope for which the survival 
probability remains asymptotically non-zero is characterized by the condition ^/a ~ 
from which we get 


1 3 

a{t) ~ 2 2 


(39) 


The term with iterated logarithm is now sub-leading, so it is asymptotically negligible. 
Thus the envelop grows according to Eq. (15). To contrast with the Brownian case, (15) 
can be called the law of the simple logarithm. 
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The trajectories of the tagged particle in single-hle diffusion can be characterized at a 
macroscopic (hydrodynamic) scale using the macroscopic fluctuation theory (MFT) (2^ . 
Indeed, the single-hle constraint allows one to express the displacement of the tagged 


particle through the density held 11,12 and thereby to apply the MFT to analyze the 


dynamical properties of the tagged particle. 

Recall that the starting point of the MFT is a hydrodynamic equation (or Langevin 
equation) for the huctuating macroscopic density prohle p{x,t): 


dtp = dr D(p)drP+ \/ ct ( p ) i ) 


(40) 


where pixfi) is a Gaussian noise with mean zero and covariance 

{ri{x, t)ri{x', t')) = 6{x — x')5{t — t'). (41) 

The angular bracket denotes ensemble average. On the macroscopic scale, all the 
microscopic details of the dynamics are encapsulated into the transport coefficients 


D{p) and cr(p), the dihusivity and mobility, respectively 30 


Consider the trajectory of the tagged particle X{f) in the time window [0,T]. The 
position of the tagged particle at any time is related to the macroscopic density prohle 
by the condition that total number of particles on either side of the tagged particle is 
invariant over time. This leads to 

pXp) POO 

/ p{x,t)dx = / [p{x,t) — p{x,0)]dx (42) 

Jo Jo 

where we have assumed that the tagged particle starts at the origin at t = 0. 

Thus the tagged particle position is a functional X(t) = Xfip] of the density 
prohle p{x,t). The statistics of the trajectories X{t) can be characterized in terms 
of the generating functional (exp[J dt X{t)X{t)]). In our previous work (n,12 we had 
X{t) = X6{t — T), so we studied the generating function (exp[AX(T)]) encoding the 
full statistics of the displacement X(T) at a single time. The general case when X{t) 
is an arbitrary function encodes the statistics over the entire time window. It is often 
convenient to study the logarithm of the average of the generating functional 

rT 


p[X,T] = log {exp 


dt X{t)X(t) 


uo 


(43) 


When X{t) = X5{t — T), the function p{X,T) = log (exp [AX (T)]) is known as the 
cumulant generating function since its expansion gives all the cumulants (X"'(T))c. 
Generally we shall call (43) the cumulant generating functional. The Legendre 


transformation of functional p[X,T] is the action A[X] characterizing the probability 
weight of each trajectory X{t) of the tagged particle, in the large T limit. 


Starting with the fluctuating hydrodynamic equation (40) we write the moment 
generating functional as a path integral 

rT 


exp 


dt X(t)X(t) 


uo 


= '^[p,p]e 


-St[p,P 


(44) 
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with action 


St[p,p]=- dtX{t)Xt[p] +F[p{x,0)] 
Jo 


+ 


dt 


dx < pdtp 


[d^p] + D{p) (4p) {da,p) \ . (45) 


This action has been derived nsing a Martin-Siggia-Rose formalism The derivation 
of (45) is essentially identical to the derivation ll |^ of the action in the special case 
when \{f) = X6{t — T), the only distinction is in the first term on the right-hand side of 
(45). Hence we jnst explain the meaning of p and F[p{x, 0)]. The field p is a Lagrange 
multiplier ensuring the validity of Eq. (40). For the deterministic initial condition the 
second term on the right-hand side of (45) is absent, so we can set F = 0; for the 
annealed setting F[p[x, 0)] is the free energy associated with the initial density profile 
p(x, 0). For the initial state with uniform average density p, this free energy can be 
written as 


F [p(a;)] = / dx 


0 V / 


(46) 


The action St[p-iP\ grows with increasing time and at large T the path integral is 
dominated by the path that corresponds to the least action. We denote this optimal 
path by (g,p) = (p, p) and perform a small variation p ^ q + 5p and p ^ p + 5p around 
this optimal path. To the leading order, the change in action ^^[p, p] corresponding to 
this variation is 


5St = 






+ 


+ 


5g(a;,0) 


5p(x,0) 


dx p(x, T) 5p(x, T) 


' —OO 

rT 


dt 


dt 


' —OO 
/•OO 


dx < —dtp 




D{q)d^a:P - X{f) 


5Xt[q] 


5q{xfi) 

dx [dtq + {(^{q)dxP) - {D{q)d^q)} 5p{x, f) 


5p{xfi) 

(47) 


where the functional derivatives are taken at the optimal path {q,p). In the above we 
have used the property that ^^[p] is a functional of only the initial profile p(x,0) and 
the profile p{x, t) at time t and has no dependence on intermediate times. 

For the action S'tIq'jP] to be minimum, the variation 6St must vanish: 5S')r[g,p] = 0. 
Ensuring this condition we arrive at the Euler-Lagrange equations and the boundary 
conditions. More precisely, since 6p{x,t) and 6p{x,t) are arbitrary, the vanishing of the 
last two integrals on the right-hand side of (47) leads to the following Euler-Lagrange 
equations 


d,q - dr (D(q)drq) = -dr {a{q)drP) 

d,p+ D(q)drrP = (drpf “ MO ' 


(48) 

(49) 
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Using Eq. (42) we rewrite Eqs. (48)-(49) as 


dtq - {D{q)d^q) = -d^ {(y{q)dxP) 

dtp + D{q)dx^p = -^cr'(g) {d^pf - B{t)<d{x -Y{t)) 


where we have used the shorthand notation 




\(t) 


q(Y(t),t) 


(60) 

(51) 

(52) 


The boundary conditions result from the vanishing of the first two integrals in (47). 
In the annealed setting both p(a;, 0) and p(x, T) are fluctuating, equivalently, 5p{x, 0) 
and 6p{x,T) are arbitrary. Then, for the minimal action, their coefficients in (47) must 
vanish, leading to the boundary condition 


p{x,0) = - 
p{x,T) = 0. 


6q{x,0) ' 6q{x,0) 


+ 


Recalling (42) and (46), we rewrite (53) as 
p(a:, 0) = 0(a:) [ dt B{t) 


2£>(r) 


a{r) 


(53) 

(54) 

( 66 ) 


For the deterministic initial condition Sp{x,0) = 0, so the first integral in (47) 
identically vanishes. The vanishing of the second integral in ( [47| ) yields ( [5^ . 

The cumulant generating functional is given by p[X, T] = —ST[q,p]- We can further 
simplify ( [45| ) using the Euler-Lagrange equations (48)-(49) to yield 

pT pT poc 


p[\,T]= dt\{t)Y{t) - F[q{x,0)] 


dt 


dx 


(T[q{xfi)] 


{d^p{x,t)) . (56) 


In deriving (56) we have taken into account the fact that the derivatives dxP and dxq 
vanish as x —)■ ±cx). 

As a self-consistency check we verify that for X{t) = XS{t — T) we recover the 
results describing the tagged particle at hnal time. 

It is instructive to compare the Euler-Lagrange equations (50)-(51) with Euler- 
Lagrange equations which appear in analysis of large deviations involving single-time 
characteristics. Most of previous studies were devoted to single-time quantities with 
problems ranging from large deviations of the stationary profile 32,33 or integrated 
current 34 in diffusive systems to the melted area in an Ising quadrant 35 and the 


survival of the target 36 . In all these examples the Euler-Lagrange equations are same. 


viz. they are identical to (50)-(51) with B{t) = 0. These equations are already very 
complicated and time-dependent solutions were found only in the simplest case of the 
Brownian single-hle system (see e.g. [^[^[^). In the present case when B{t) ^ 0 we 
do not know how to solve (g-iig even in the Brownian case. 

It is important to analyze the scaling of p[X,T] with T. For this, we rescale time 
r = t/T and position i = x/ a/T. Substituting this in the formal solution (56) and using 
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the Euler Lagrange equations (50)-(51) along with their boundary conditions we hnd 
the scaling relation 

'X{t) 




T 


-,T 


= Vf fx[X{T),l] 


(57) 


which leads to the announced sub-diffusive scaling (18) of the cumulants. 


5. Two-time Correlation 

The formalism of the previous section applies to any single-file system. It is impossible, 


however, to solve the Euler-Lagrange equations (50)-(51) in the general situation when 


the transport coefficients D{p) and a{p) are arbitrary. The framework presented in 
Sec. 1^ is still helpful to derive some general results. In this section we derive the two- 
time correlation function, 0 and ([^, for an arbitrary single-file system. We employ a 
perturbation approach which has already proved useful in a number of problems. 
In the present case, we write 

X{t) = e h{t) (58) 

and treat e as a small parameter. When e = 0 the solution is q{x, t) = p and p{x, t) = 0. 
Generally we write a series expansion 

q{x, f) = p + e qi{x, t) + q 2 {x, t) H- (59) 

p{x,t) = epi{x,t) X-P 2 {x,t)-\ -. (60) 

The cumulant generating function can also be written as a series 

p[eh,T] = e^P 2 [h,T] -. (61) 

The terms in odd powers of e vanish because the microscopic dynamics is unbiased. 


From the definition of the cumulant generating functional (43) we have 

Pn[h,T] = ^[ dti... I dtn firi) 


n\ 


(62) 


'0 


'0 


where kn{ti, • • •, t^) is the n-time cumulant. The cumulants are related to the multi-time 
correlations. In the second order, we have k 2 (ti,t 2 ) = (^(^i)-^(G))- 


Combining the series expansion (59)-(60) and the Euler-Lagrange equations (50) 


(51) we find 


dtqi - D{p)d^^qi = - a{p)d^^pi 

h{t) 


dtpi D{p)da,xPi = 


P 


-0(a:) 


in the first order. In deriving Eqs. (63)-(64) we have taken into account that 


B[t) = PXXl + o[r). 

p 


(63) 

(64) 

(65) 


This can be seen from the definition of B{f) in (52) 
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We write a similar series expansion for the fnnction Y(t) 
Y (t) = eYfif) + e^Y2{t) + ■ ■ ■. 


Using (42), we establish the vanishing of the zeroth order term and find 


1 /■°° 

^i(^) = - / dx[qi{xfi) - qfix,^)] 
P Jo 


Snbstitnting all these series expansions into (56) we get 

qi2[h,T]= f dt h{f)Yi{f) - F 2 - / dx{d^pif 

Jo ^ Jo J -00 

with F 2 = 0 for the deterministic initial condition; in the annealed case 

= [ dx{qi{x,Qi)f 

^\P) J — 00 


( 66 ) 

(67) 

( 68 ) 

(69) 


follows from Eq. (46). We now determine {X{ti)X{t 2 )) by solving Eqs. (63) 


5.1. Deterministic initial condition 
The bonndary conditions are 

Pi(a;, T) = 0 and gi(a;,0)=0. 


The solution of (64) subject to pi{x,T) = 0 reads 

rT roo 


Pi{x,t) = p / dr dz h{T) Q{z) gD{z,T\x,t) 
J t J —00 


where 


gD{z,T\xfi) = 


: exp 


[z — X] 


,jATfJD{p){T^Ff^ L 4T)(p)(r-t)_ 
for 0 < f < r. Performing the Gaussian integral over z we get 


Pi{xfi) = — / dr h(r)erfc I 




2pi. . '''-\2pD(p)(T-t) 


We shall also need dxPi, given by 


%Pi(x,t) = - / dr h{T) go{0,T\x,t) 

P Jt 


(70) 


(71) 


(72) 


(73) 


(74) 


It is convenient to write qi(x.t) — Plugging this into (63) and solving 

the resulting equation we obtain 

pt poo 

fi{x,t) = a{p) dr dz gD{x,t\z,T) d^pi{z,T). (75) 


We now calculate P 2 [h]det for the deterministic initial condition, F 2 = 0. Eq. (68) 
becomes 


hit) aio) 7“ 

h 2 [h,T]det= / dt - / dx qi{x,t) -— dt dx {d^Pif 

Jo P Jo ^ Jo J -00 

= h-h (76) 
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where we have used Eq. (67) and q'i(x,0) = 0. 


To compute we write again qi = integrate by parts, and take into account 

that vanishes at a: —>■ oo: 


1 

Ji = - / dt hit) 
P Jo 


(77) 


Taking from Eq. (75) and d^pfiz^r) from Eq. (74), we re-write as 


0-( q\ p^ 

h = —^ / dt hit) dr dti h(ti) / dz gD{0,ti\z,T)gD{0,t\z,T). (78) 

P Jo Jo Jt J—oo 

Computing the Gaussian integral over z, we arrive at 

I it h(t) I dr I dt, Kt0^j^D= (79) 

where we used the shorthand notation (|^. The integration over r variable is performed 
using identities like 

P dr 


= y/ti -\-t — \/ti — t when ti > t. 


Jo y/ti + ^ — 2 t 

This leads to 

rT rT 

fi = V{p) / dti / dt2 h(ti)h{t2) y/ti -\-t2- \/\ti - t2 
Jo Jo 

where we denote t 2 = t. 


(80) 


(81) 


To compute I 2 , we take dxPi{x,t) from Eq. (74) and obtain 


0 -(p\ p'^ p'^ ^00 

h = ^ r, / dt dr dti h{T)h{ti) / dx (7£)(0, fi|a;, f)5fD(0, r|x, t). ( 82 ) 

• 0 Jt J t J —00 


2 p 2 


Changing the order of integration, dt dr = Jq dr dt, we find J 2 = |/i once we 


recall Eq. (78). Therefore 


/X2[h,T]det = / dti [ dt2 h{ti)h{t2) Vh + t2 - \/|fi - t2\ (83) 


Comparing this with Eq. (62) for n = 2 we arrive at the announced result for the 
two-time correlation 


{X{ti)X{t2))det — T^ip) \/tl + ^2 — \/\tl — 


(84) 


5.2. Annealed setting 

The boundary conditions for pi and qi are now given by 

<^{p) 


Pi{x,T) = 0, and gi(x, 0) = 


2D{p) 


Pi(x,0) — Q(x) [ dt^^^^ 

Jo P 


Using (67)-(69) we get 


1 


/i2[h,T]ann = - / dt h{f) / dx [qfix, t) - qfix, D)] 

P Jo Jo 

f*T POO 

dt 


D{p) 


(^ip) J-c 


dx {qi{x,0)f - 


dx {dxPi{x,t)y 


(85) 


( 86 ) 
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In calculations, it is convenient to focus on fi 2 [h,T]s.nn — t^ 2 [h,T]det- For both the 
deterministic and the annealed settings, the governing equation for pfixfi) and the 
boundary condition on pfix^T) are identical, so the solutions coincide, and are given 
by Eq. (73). The governing equation (63) for qfixfi) is also the same in both cases, but 
the boundary conditions are different. Since Eq. (63) is linear in qi{x,f), we write the 
solution as a sum 


qi{x,t) = qfixfi) + qh{x,t) 

where qfixfi) is the solution of the inhomogeneous diffusion equation 
dtqi - D{p)d^^qi = -(y{p)d^^pi with qfix, 0) = 0 
and qh{x,t) is the solution of the homogeneous equation 
dtQh dD (^p^OxxQh 0 

subject to the boundary condition 


Qh{x,0) = 


2D{p) 


Pi{x, 0) — 0(x) 


dt 


i‘(t) 

p 


(87) 

( 88 ) 

(89) 

(90) 


Since qfixfi) and qi{x,t) in the deterministic case satisfy the same boundary-value 
problem, they are identical. Further, using ( |7^ and (86) we conclude that 

r*T /*oo 




1 

P 
D{p) 


P 2 [h,T]det = - / dt hit) / dx[qhix,t) - qhix,0)] 
P Jo Jo 


dx iqhix,0))'^ . 


(91) 


^{p) J-c 

Thus the difference depends only on qh{x,t). Further simplihcation comes from the 
identity (proved in Appendix A): 


P Jo 


dt hit) / dx [qhix, t) - qhix, 0)] = 


2 Dip) r 

(^ip) J-c 


dxiqhix,0))‘^ . (92) 


Thus, 

f^2[KT], 


/^2 ^]det 


Djp) 

(^{p) J- 

^(p) 

^Dip) 


dx iqhix,0)y 


OO 

OO 


dx 


Piix, 0) — ©(x) 


dt 


hjt) 

P 


(93) 


where in the last step we have used the boundary condition 

J i; 

^(p) 


We extract pi(a;, 0) 


from (73) and after straightforward calculations recast (93) into 
P2[h,T]ann - P2[h,T]det = / dti dta h(ti)h(t2) 


SpWip) 


[ dx erfc (^ erfcf , J] 

Jo \o/mwp \vwphj 


Finally, using the identity 38 

poo 

/ dy erfc 


(w) “K 


y 

a/^ 


— -J= \\/D + \ft2 — Vd + ^2 ] 

Jtt 


(94) 


(95) 
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we perform the integral over x and reduce Eq. (94) to 

rT pT 


f'2[h,T]ann — f'2[h,T]det — —I dti I dt2 h{ti)h{t2) [y/tf + \/t2 — + ^2 ] 

^ Jo Jo 

which in conjunction with Eq. (83) lead to the announced result 


{X(ti)X{t2))ann — dJip) VD + ~ \/|^l ~ ^2 


(96) 


6. Summary 


We have investigated dynamical properties of a tagged particle in single-file diffusion. We 
computed the two-time correlation function for the displacement of the tagged particle. 
The results ([^-([^ are simple, but the derivations are involved since they rely on the 
macroscopic fluctuation theory (MET). The advantage is the generality: our predictions 
(§-(§ apply to an arbitrary single-file system, all the specificity of a concrete system 
is encapsulated in an amplitude D{p) which is a simple combination of the fundamental 
transport coefficients, D[p) and cr(p), characterizing the concrete single-hle system. For 
a single-file system of Brownian particles with hardcore repulsion, our results (|^-(|^ 
are confirmed by an independent microscopic calculation j^. 

Using the two-time correlation function, we calculated two non-trivial moments 
of the fraction of time which the tagged particle spends to the right of its starting 
point. The major challenge is to compute all such moments, equivalently to find the full 
probability distribution. For Brownian motion, this probability distribution is governed 
by the arcsine law. We also discussed a single-file analog of the law of the iterated 
logarithm describing the minimal asymptotic envelope of the Brownian motion. This 
analog, the law of the ‘simple’ logarithm, asserts that the ratio of maximal fluctuations 
of the tagged particle to typical fluctuations grows as v^lnt. 

Using the MET we devised a formalism giving the cumulant generating functional. 
The cumulant generating function containing the full single-time statistics of the 
displacement of the tagged particle is a very special outcome of the formalism. In 
principle, the formalism gives the full statistics over the entire history, so all multiple¬ 
time correlation functions are encapsulated in the cumulant generating functional. As 
it often happens, the chief problem is that the governing Euler-Lagrange equations 
which one needs to solve to extract concrete results are extremely complicated. In 
other problems relying on the MET, essentially the only solvable case corresponds to 
Brownian particles with hard exclusion constraint. In the present case even this seems 
to be very non-trivial, so far the solution for the generating function is known only in 
the single-time situation. 

Finally, we mention that the Brownian motion has a number of remarkable 


properties, such as beautiful laws 19 21 governing its maximum. Single-hle analogs of 


these properties are unexplored and fascinating open problems. 


We are indebted to Davide Gabrielli for asking questions that have triggered this work. 
TS thanks Bernard Derrida for very useful discussions on the single-hle problem. This 
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Appendix A. Derivation of the identity (92) 


Denote by L the left-hand side of Eq. (92). We start by rewriting the expression as 

'hit) 


L = 


dx / dt 


P 


-0(a:) 


(Ihixfi) 


r 

dx / dt 


'0 


hjt) 
L P 


0(x) 


g/*(x,0). (A.l) 


Using the boundary condition (90) and performing a rearrangement of the terms we get 
2D(p) 


with 


L = 


U = 


(^ip) J-O 


dx[qh{x,0)] -V 


dx 


pi(x,0)gfe(x,0) - p 


-1 


dt h{t)Q{x)qh{x, t) 


'0 


(A.2) 


(A.3) 


To establish Eq. (92) it suffices to show that U = 0. Using Eq. (64) we re-write V as 

/ oo poo pT 

dxpi{x,0)qh{x,0) + dx dt qh{x,t) [dtPi{x,t)-P D{p)d^^pi{x,t)]. 

■OO J — oo J 0 

Integrating dtqh{x,t)dtPi{x,t) by parts and recalling that pi{x,T) = 0 we get 


U = 


dx dt[Dip)qh{x,t)d^^pi{x,t)-pi{x,t)dtqh{x,t)]. (A.4) 


Using (89) we reduce (A.4) to 


U = Dip) 


dt 


dx dx [qh{x, t)dxPiix, t) - pfix, f)dxqh{x, t)] (A.5) 


which clearly vanishes. Thus U = 0 as we claimed. This completes the proof of (92). 


References 

[1] Richards P M 1977 Phys. Rev. B 16 1393 

[2] Karger J and Ruthven D 1992 Diffusion in zeolites and other microporous solids (Wiley) 

[3] Kukla V, Kornatowski J, Demuth D, Girnus I, Pfeifer H, Rees L V C, Schunk S, Unger K K and 

Karger J 1996 Science 272 702 

[4] Li G, Berg O G and Elf J 2009 Nature Physics 5 294 

[5] Das A, Jayanthi S, Deepak H S M V, Krishna V R, Kumar A, Dasgupta G and Sood A K 2010 

ACS Nano 4 1687 

[6] Siems U, Kreuter G, Erbe A, Schwierz N, Sengupta S, Leiderer P and Nielaba P 2012 Sci. Rep. 2 

1015 

[7] Harris T E 1965 J. Appl. Probab. 2 323 

[8] Arratia R 1983 Ann. Probab. 11 362 

[9] Rodenbeck C, Karger J and Hahn K 1998 Phys. Rev. E 57 4382 

[10] Leibovich N and Barkai E 2013 Phys. Rev. E 88 032107 

[11] Krapivsky P L, Mallick K and Sadhu T 2014 Phys. Rev. Lett. 113 078101 



















Dynamical properties of single-file diffusion 


20 


[12] Krapivsky P L, Mallick K and Sadhu T 2015 J. Stat. Phys. 160 xxx 

[13] Sethuraman S and Varadhan S R S 2013 Ann. Prob. 41 1461 

[14] Hegde C, Sabhapandit S and Dhar A 2014 Phys. Rev. Lett. 113 120601 

[15] Wei Q H, Bechinger C and Leiderer P 2000 Science 287 625 

[16] Lutz C, Kollmann M and Bechinger C 2004 Phys. Rev. Lett. 93 026001 

[17] Lin B, Meron M, Cui B, Rice S A and Diamant H 2005 Phys. Rev. Lett. 94 216001 

[18] Krug J, Kallabis H, Majumdar S N, Cornell S J, Bray A J and Sire C 1997 Phys. Rev. E 56 2702 

[19] Ito K and McKean H P 1965 Diffusion processes and their sample paths Die Grundlehren der 

mathematischen Wissenschaften (Berlin: Springer) 

[20] Feller W 1968 An introduction to probability theory, Vol. 1 3rd ed (New York, NY: Wiley) 

[21] Morters P and Peres Y 2010 Brownian Motion (New York: Cambridge University Press) 

[22] Kac M 1949 Transactions of the American Mathematical Society 65 1 

[23] Majumdar S 2005 Current Science 89 2076 

[24] Bertini L, De Sole A, Gabrielli D, Jona-Lasinio G and Landim C 2007 J. Stat. Mech. P07014 

[25] Derrida B 2007 J. Stat. Mech. P07023 

[26] Bertini L, De Sole A, Gabrielli D, Jona-Lasinio G and Landim G 2014 ArXiv e-prints 1404.6466 

[27] Abramowitz M and Stegun I A 1964 Handbook of mathematical functions with formulas, graphs, 

and mathematical tables vol 55 (Washington, D.C.: U.S. government printing office) 

[28] Bacon R H 1963 Ann. Math. Statist. 34 191 

[29] Benichou O and Desbois J 2015 J. Stat. Mech. P03001 

[30] Bodineau T and Derrida B 2004 Phys. Rev. Lett. 92 180601 

[31] Martin P, Siggia E and Rose H 1973 Phys. Rev. A 8 423 

[32] Tailleur J, Kurchan J and Lecomte V 2007 Phys. Rev. Lett. 99 150602 

[33] Bertini L, De Sole A, Gabrielli D, Jona-Lasinio G and Landim C 2001 Phys. Rev. Lett. 87 040601 

[34] Derrida B and Gerschenfeld A 2009 J. Stat. Phys. 137 978 

[35] Krapivsky P L, Mallick K and Sadhu T 2015 J. Phys. A 48 015005 

[36] Meerson B, Vilenkin A and Krapivsky P L 2014 Phys. Rev. E 90 022120 

[37] Krapivsky P L and Meerson B 2012 Phys. Rev. E 86 031106 

[38] Prudnikov A P, Brychkov Yu A and Marichev O I 1986 Integrals and series: special functions vol 2 

(New York: Gordon & Breach Science Publishers) 

[39] Sadhu T and Derrida B 2015 J. Stat. Mech. xxx xxx 



